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G> ■ Abstract 



In the leading 0{ag) order of the perturbative QCD, the hadronic production 
D ' cross-sections of S- and P-wave states of Be meson are calculated. The results for 

the S-wave levels are compared with the values given in other papers as well as in 
, the model of the b quark fragmenatation into Be- In the given order, the cross- 

I sections of the hadronic production of the P-wave states are calculated for the first 

' time. Their contribution into the Be meson production is less than 10%. There 

. is a strong difference between the predictions of the fragmentation model and the 

I exact perturbative calculations. These difeerences are discussed in details for the 

Cf^ ■ differential distributions over variuos kinematical quantities. 

(N : 
o . 

^ ■ Introduction 

CL|| During recent two years an essential progress has been reached in the understanding of 

p • the production mechanisms for the heavy quarkonium, composed of two heavy quarks 

^ ■ with different flavours. 

The 6c-quarkonium production process has the most simple consideration in e^e~- 

! annihilation, where in the limit of high energies {M'^/s -C 1), the differential cross-section 

^ ' has the evident factorizable form 

■ - - ' 0,(7 ]-[ / \ 7-1 / \ / \ 

= (^bb{s) ■ Di_,h{z), (1) 

where z = 2Eb^/ a/s, and Di^jj{z) is interpreted as the b — H-\-X fragmentation function 
with H being the 6c-quarkonium [1-3]. 

Recently, the functions of fragmentation into the 5*-, P-, and Z^-wave states for the 
6c-quarkonium production in e^e~ collisions have been found [1-4]. Everywhere below 
mentioning the fragmentation model we will imply expressions like (0). 

The description is more complicated for 77 and hadron-hadron collisions. As for these 
interactions, in addition to the mechanism of the b quark production with the following 
fragmentation into the 6c- quarkonium, there is a new type of diagrams, corresponding to 
the photon (gluon) dissociation into the pair of heavy quarks with the following recombi- 
nation of the quarks into the 6c-quarkonium. 

In the photon-photon collisions, one can separate three gauge invariant subgroups 
(6-|-6-|-8) within complete set of 20 diagrams in the leading Born approximation, so that 
subgroups allow one to have the self-consistent interpretations: the b quark fragmentation. 
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c quark one, and recombination. The contribution of the first mentioned type is quite 
rehably described by expression (|ID. For the c quark fragmentation diagrams one observes 
a strong deviation from the picture of the fragmentation model. The recombination 
contribution dominates in the complete region of kinematical variables. As was shown in 
[5], the given conclusion is valid for the 5*- wave state production as well as the P-wave 
one [6]. 

The process of hadronic production of the 6c-quarkonium is more difficult for the 
analysis. In the leading order of the perturbation theory it is described by 36 diagrams 
in the fourth order over a^. At present, there are results of calculations by several groups 
giving controversial numerical values and conclusions for the 6c-quarkonium production 
cross-sections [7-10]. Moreover, it was offered to simplify the consideration of the heavy 
quarkonium production by its reduction to the straightforward fragmentation of b quark 
with the usage of the fragmentation functions, derived for e"'"e~-annihilation case [11]. 

In this paper we give the comparative analysis of the results by different papers devoted 
to the production of the S'-wave states. We prove that the straightforward fragmentation 
model does not properly describe the 6c-quarkonium production. Next, we present the 
results of calculations for the P-wave level production over the complete set of the leading 
order diagrams. 



1 Calculation technique 

The A^"^^^ amplitude of the Be meson production can be expressed through the amplitude 
of four free quarks production T^^^ {pi, fc(q)) and the orbital wave function of the meson, 
\I'-^'^(q), in the meson rest frame as 



where J and jz are the total spin of the meson and its projection on z axis in the Be rest 
frame, correspondingly; L and Iz are the orbital momentum and its projection; S and Sz 
are the sum of quark spins and its projection; C^-'/^ are the Clebsh-Gordan coefficients; pi 
are four-momenta of Be, b and c, q is the three-momentum of b quark in the Be meson 
rest frame; k{q) is the four-momentum, obtained from the four-momentum (0, q) by the 
Lorentz transformation from the Be rest frame to the system, where the calculation of 
T^^'^ipi, k[c\)) is performed. Then, the four-momenta of b and c quarks, composing the Be 
meson, will be determined by the following formulae with the accuracy up to |qp terms 

Pc = ^PB.-Mq), ^' 

where m?, and me are the quark masses, M = mf, -|- me, and Pb^ is the Be momentum. 
Let us note that for the P-wave states it is enough to take into account only terms, linear 
over q in eq.(^, and q = in the S'-wave production. 

The product of spinors viUc, corresponding to the b and c quarks in the T^'^^i^pi, k{q)) 
amplitude of eq.(|D, should be substituted by the projection operator 

. X / / TT^Bc + k-mb\ I ^Pb, -k + me. 
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where F = 7^ for = 0, or F = i*{PB^, Sz) for S = 1, where s{Pb^-, Sz) is the polarization 
vector for the spin-triplet state. 

For the sake of convenience, one can express the V{T) operator through the spinors of 
the following form 



u'^iPc -k,±) = [l-^J u,{p„ ±), 

where Vh{phi ±) and Uc{pc, ±) are the spinors with the given projection of quark spin on 
z axis in the meson rest frame. Note, that the spinors in eq.(^ satisfy the Dirac 
equation for the antiquark with the momentum Pb + k and mass nib or for the quark with 
the momentum pc — k and mass rric up to the linear order over k ( i.e. over q, too), 
corresp ondingly. 

One can easily show that the following equalities take place 



2M 1 



2m,2m, Tl^^^lPfe + k, +X(Pc - k, +) - v',,{pb + k, - k, -)} 



2M 



2m,2mybiPb + k,+HiPc-k,-) 

= V{e*{P,-l)) + 0{e), 



(6) 



2M 1 



2m,2m, V2^'"biPb + k, +)u'^{Pc -k,+)+ v[{pb + fc, -)u'^{Pc -k,-)] 

= v{i*{p,o)) + o{e), 



^^v',{pb + k,-)u',{p,-k,+) = 
= V{e*{P,+l)) + 0{e). 

In the Be rest frame, the polarization vectors of the spin-triplet state have the form 

£^/(-l) = ^(0,l,-z,0), 
£'-/(0) = (0,0,0,1), 
^'■^(+1) = -71(0,1,^,0). 

In calculations the Dirac representation of 7-matrices is used and the following explicit 
form of the spinors is applied 



(7) 



u{p, +] 



I E + m \ 




Pz 



u{p, -) 



\ Px+ iPy / 





E + m 
Px - iPy 
-Pz 



v{p, +] 



/ Pz 
Px + iPy 



\ E + m J 



v{p, 



( Px- iPy \ 

Pz 



V E + m J 



The S'-wave production amplitude can be written down as 

T 



A' 



2M 



2mh2mr V 47r 



T'^' {p^,k{q = 



(9) 
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where Rs{0) is the radial wave function at the origin, so that 



Rs{0) 



2 JSc! 



and the Jb^ value is related with the leptonic constants of pseudoscalar and vector Be 
states 



where J^(x) and J^^{x) are the vector and axial- vector currents of the constituent quarks. 
Then the account for hard gluon corrections in the first order over [12] results in 



/ 
/ 



fv 
fp 



/ 1712 — mi ^ ^2 8 
1 In 

TT \m2 + mi mi 3 



1 



a 



H 



m2 - mi m2 

m 2 

m2 + mi mi 



(10) 
(11) 



where mi^2 are the masses of quarks, composing the quarkonium. For the vector currents 
of quarks with equal masses, the BLM procedure of the scale fixing in the "running" 
coupling constant of QCD [13] gives (see paper by M.B.Voloshin in ref. [12]) 



a. 



H 



^MS. -11/12 . 



The estimates of the fs^ value within the potential models have the essential uncertainty, 
/bc = 500 ± 100 MeV [14]. The QCD sum rule estimate of the fs^ value for the pseu- 
doscalar state gives fs^ = 385 ± 25 MeV [15], which is in a good agreement with the 
evaluation in the framework of recent lattice computations [16], where /^^ = 395(2) MeV 
with the error bar, giving the statistical uncertainty only. The /^^ estimate strongly de- 
pends on the scale choice, which is not yet calculated in the BLM procedure. So, we 
use Jb, = 570 MeV. 

For the P-wave states in eq.(|^), the T^'^^ {Pi,k{q)) amplitude can be expanded into 
the Taylor series up to the terms linear over q. Then one gets 



;i2i 



where -Rp(O) is the first derivative of the radial wave function at the origin, and has 
the following form 

^° = W.^ (13) 

where ^ are the differential operators acting on T^'^' (pi, k(q)) as the function 

of q = {qx, qy, qz) at q = 0. 
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As all considered matrix elements are calculated in the system distinct from the Be 
rest frame, the four-momentum A;(q) has been calculated by the following formulae 



7,0 _ v q 

k = q+(7T=^-l)^v, ^^"^ 

where v is the Be velocity in the system, where the calculations are performed. The 
matrix element T"^*^ {pi, k{q)) is computed, so that the four-momenta of b and c quarks 
are determined by eq.(H), taking into account eq.(|T4|). 

The first derivatives in eq.(|13D are substituted by the following approximations 

% '''=°^ A ' ^^^^ 

where A is some small value, and q-^ have the following form 

q- = (A,0,0), 

q^=(0,A,0), (16) 
q^ = (0,0,A). 

With the chosen values of quark masses and interaction energies, the increment value 
A = 10^^ GeV has provided the stability of 4-5 meaning digits in the squared matrix 
elements summed over jz for all P-wave states with the given value of J and S, when one 
has performed the Lorentz transformations along the beam axis or the rotation around 
the same axis. 

One has to note that because of such transformations, the new vectors k{q^) do not 
correspond to the transformed old vectors. Therefore, the applied test is not only a check 
of the correct typing of the T^'^' {pi,k{q)) amplitude, but it is also the check of correct 
choice of the phases in eq.(lT^). 

The matrix element A^'^^'' squared, which is calculated by the method described above, 
must be summed over as well as the spin states of free b and c quarks. It also must be 
averaged over spin projections of initial particles. 

The phase space integration has been made by the Monte Carlo method of RAMBO 
program [11]. 



2 S'-wave states 

Let us consider the calculation results for the (7(7-production of the S'-wave 6c-quarkonium 
levels {Be and B*) in comparison with the values given in other papers. The cross-section 
of the process under consideration is proportional to the /^^ squared as well as to the 
fourth power of a^. To compare the results of different papers, we rescale all the numbers 
to the values, determined by the same set of = 0.2 and /^^ = 570 MeV, which are used 
in the given calculations. We also fix the mass values of = 4.8 GeV and rric = 1.5 
GeV. 

The cross-sections of the gluonic production of 5^*^ versus the total energy of the gg 
collisions are presented in Fig. 1. In addition to the results of the complete numerical 
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calculations of the 0{a'^) contribution, we give also the values, obtained in the fragmen- 
tation model and, hence, calculated as the product of the (y^^f-production cross-section for 
the 66-pair and the probability of the b —>■ B^*^ fragmentation, calculated at the same set 
of parameters. 

As one can see in this figure, there is a good agreement of our previous calculations 
[7] with the results of refs.[8,10]. The values, given in two other papers of ref.[9], are 
approximately three times greater than our results. In this paper we have recalculated 
the cross-sections in the axial gauge, used in ref. [9], in contrast to the covariant Feynman 
gauge applied in our previous consideration [7]. The corresponding numerical results 
have not changed after the replacement of the gauge. This fact points to errors in the 
computations, performed in ref. [9]. 

In the studied region of energies, the ratio of the Be and B* yields, R = (Jb*/(JBc — 3, 
strongly deviates from i? ~ 1.4 predicted by the fragmentation model. One can also see 
in Fig. 1, that in the region of the applicability determined by the condition M^/s -C 1, 
the fragmentation model gives the total cross-section, which is essentially smaller than 
the exact result evaluated over the complete set of diagrams. 

As one can see in Fig. 2, the total cross-section of the S^*^ production is basically 
accumulated in the region of the transverse momenta close to pt ~ Mb^- One could 
expect, that the fragmentation mechanism begins to dominate at large pt- Indeed, the 
PT-distribution shown in Fig. 2 for Be and B* at = 100 GeV, points to the fact, that 
there is a quite good coincidence of the distribution tails obtained from exact perturbative 
calculations and in the fragmentation model. Numerically, this agreement takes place at 
Pt > 40 GeV for Be- Thus, one can draw the certain conclusion, that the fragmentation 
contribution is not dominant in the (yff/- collisions. It also does not describe the distribution 
over the variable defined as 2; = 2\PbJ/\/§- This fact is evident in Fig. 3. At z values 
close to unit, the fragmentation model overestimates the exact perturbative result. 

To get the cross-section and px spectra in hadron-hadron collisions, one should con- 
vert the partonic gg cross-section with the distribution functions of gluons in the initial 
hadrons. We use the parametrization of ref. [17] for the parton distributions at the fixed 
virtuality scale Q — 2mi, ~ 10 GeV. In the framework of the given approximation, we fix 
also the coupling constant value ag — 0.2. 

The convolution result is presented in Fig. 4 for the energy of the FNAL Tevatron 
{^/s = 1.8 TcV). The histograms correspond to the exact perturbative results, whereas 
the curves are calculated according to the following formula 

H{pT)x) = Y. dxidx2dzfi/p{xi, n)fj/p{x2, n) x 

b{pT/z) + x)xDi^H{z,IJ'), (17) 

where D(z, ji) is the function of the 6 — > if fragmentation with H — B^ -B*..., da/dpT is 
the differential cross-section of the h quark production, and fi/A{x,ii) is the distribution 
of the i-kind parton in the A-hadron. 

As one can see in the figure, the curves certainly deviate from the histograms of the 
exact perturbative calculations in the all region of px- At small px, the fragmentation 
model gives the overestimation, whereas at large px, contrary, it underestimates the exact 



X- 



da 
dpT 
da 
dpT 
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result. It is significant to note the fact, that the perturbative calculations and the model 
give the different values for the ratio of the differential cross-sections of B* and B^. One 
finds the yield ratio i? ~ 2 3, whereas the model gives R ~ 1.3 -r- 1.5. 

The qualitative agreement between the model and perturbative calculations for the 
da/dpT distribution of Be meson has allowed the authors of ref.[8] to conclude on the 
satisfactory description of the exact perturbative 0(ag)-contribution by the fragmentation 
model. The misleading is related to the fact, that those values of the gfgf-collision energy 
are included in the integration region of expression (|1^, where the approximation of 
the fragmentation does not work or it is not strictly defined, i.e., where the M^/s -C 1 
condition is evidently not valid. As one can see in Fig. 1 at these low energies, the cross- 
section in the fragmentation model and that of the perturbative result have the different 
dependency on the energy. The latter fact is because of the use of the two-particle phase 
space in the fragmentation approximation of the heavy quark production instead of the 
real three-particle phase space {Be, b and c are in the final state). 

To avoid the uncertainties related with the low energy of ^f^^-collisions, we have started 
the integration over s in expression from \/§ > 60 GeV. As one can see in Fig. 5, the 
given cut drastically changes the relation between the fragmentation model contribution 
and the exact perturbative result. It indicates that the fragmentation model can not 
provide an adequate approximation for the correct description of the hadronic Be meson 
production. 

One has to note, that the given conclusion does not depend on the definite choice of 
the structure functions for the initial hadrons as well as on a special modification of the 
fragmentation expression. 

Thus, the conclusion drawn in ref. [8] about the dominance of the fragmentation mech- 
anism in the hadronic production of B^ mesons, is incorrect. 



3 Production of P-wave states 

As it was mentioned above, to get the cross-section for the production of the P-wave 
quarkonium levels it is necessary to calculate the first derivative of the quark production 
matrix element over q, the relative momentum of the quarks inside the Be meson. This 
procedure of calculations has been tested by comparison of the numerical computation of 
the fragmentation functions for the b quark into the P-wave states with known analytical 
expressions for the corresponding functions obtained in ref. [4]. We have found a good 
agreement between the ^Pi, ^Pq, ^Pi, and ^P2-level cross-sections calculated numerically, 
and those obtained from the analytical expressions for the same set of parameters, rrih, me, 
I P^ (0)1 2 and as. The total cross-sections of the gluonic production of P-wave levels are 
presented versus the total energy of the collisions in Tab. 1. The dependence of the cross- 
section, summed over the spin states of the P-levels, on the energy of the interacting gluons 
is shown in Fig. 6. This dependence can be approximately described by the following 
expression 

= 25, . fi - f ?i!!!±±!nA] r\( n^nip!^] " pb, (18) 



The contribution, determined by the fragmentation model, is also shown in the same 
figure, where it is found as the corresponding total gg — bb cross-section multiplied by 
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the integral probability of the b Bc{L = 1) fragmentation, W = 5.34 ■ 10~^. The 
c-quark fragmentation is suppressed by the order of magnitude in the hadron collisions, 
and we will neglect it below. 

As one can see in Fig. 6, the total cross-section of the P-wave level production is 
much greater than the value predicted by the fragmentation model, which has a strict 
meaning in the region of high energies, and, thereby, the contribution of the recombi- 
nation diagrams dominates. The fragmentation contribution also does not describe the 
da/dz distributions, which are the differential cross-sections integrated over the transverse 
momentum of P-wave levels. These distributions are shown in Fig. 7. It is interesting 
to notice that in contrast to the 77-collisions, the result of the fragmentation model es- 
sentially overestimates the exact perturbative values in the gluonic collisions. This fact, 
related with the interference of different contributions, will be discussed elsewhere. 

The distributions over the transverse momentum of the produced P-wave states are 
presented in Fig. 8 in comparison with the results of the fragmentation model. One can 
see, that the fragmentation is valid only on the tails. This fact is quite expected, since 
the same picture has been observed in the production of S'-levels. 

Thus, our calculations of the P-wave level production cross-sections show that the 
total hadronic cross-sections are about 10% of the corresponding cross-sections of the 
S'-wave levels, and the recombination mechanism dominates for the both S- and P-levels. 

The differential distributions for the P-wave states calculated by convoluting the par- 
tonic cross-sections with the gluon distributions, as it was described in previous Section, 
are shown in Fig. 9 for the FNAL Tevatron energy, ^/s = 1.8 TeV. 

As it comes for the S'-wave states, the calculation over formula (|l]) of the fragmenta- 
tion model gives a rather approximate qualitative value for the P-level cross-section. The 
reason of such "description" is generally related to the uncertainty of the fragmentation 
approach in the region of energies close to the kinematical threshold of the reaction, where 
the fragmentation model overestimates the cross-section. This incorrect contribution com- 
pensates the underestimation in the region of large energies and small or moderate pt, 
where the recombination dominates. This results in the approximate description of the 
exact perturbative distribution by the fragmentation formula at pt ~ 15 GeV (see Fig. 9). 
As it was shown in previous Section, the s cut clarifies the problem, and one evidently 
conclude, that the exact perturbative calculations and the fragmentation model strongly 
differ in the predictions. 

Conclusion 

Let us analyse the results of the numerical perturbative calculations performed for the 
hadronic production cross-sections of the S- and P-wave levels of the 6c-quarkonium in 
the leading 0{aj) order of the perturbative QCD. 

There are two scales of virtualities in the problem. The first one is of the order of the 
heavy quark masses, and it appears in the heavy quark production. The second scale is 
determined by the characteristic relative momentum of the quarks inside the 6c-meson. 
Since the latter is much less than the masses of the produced quarks, the process of the 
Be production can be separated into two steps: 

1) the heavy quark production, described in the perturbative QCD, and 
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2) the forming of the bound state, described by the quarkonium wave function at the 
origin. 

The results, hence, Uncarly depend on a^, /^^ and |i?p(0)p. Moreover, there are the 
additional parameters defining the quark masses, and rric- The bulk of the uncertainty 
of the results is connected with the q;^, /s^, and |i?p(0)| values. For instance, the use of 
the running coupling constant as{s) instead of Q;s(4m^) as well as the use of a smaller 
value of fs^ can decrease the cross-section by order of magnitude. 

From our point of view, the use of as{s) is not correct, since the analysis shows that 
typical virtualities in the production of four heavy quarks are less than s at large s. This 
is why we have fixed the value of ag — 0.2. The comparison with the results of other 
papers shows, that we agree with the calculations of the yS-level production in ref.[8], 
but, in contrast to [8], we compare the yields of pseudoscalar and vector states of B^. in 
the hadronic production to emphasize the invalidity of the fragmentation model to the 
problem under consideration. The results of ref . [9] disagree with the values given in our 
paper for both total cross-section of the pseudoscalar state and the dependence on the 
total energy. The production of the vector state was not considered in ref. [9]. So, our 
computations, performed in covariant and axial gauges, point to errors in papers of ref. [9]. 

The exact perturbative cross-sections for the P-wave states are found for the first time. 
As well as for the hadronic production of S'-levels, one observes the strong discrepancy 
in the values of the relative yields of ^Pi, ^Pq, ^Pi and ^P2-levels in comparison with the 
production in e'''e~-annihilation. This fact points to the different dominating mechanisms. 
The fragmentation of the b quark into B^*^ is the basic mechanism in e"'"e~ -annihilation, 
whereas the recombination of heavy quarks dominates in the hadronic production of B^*^ . 
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Figure captions 



Fig. 1. The total cross-section of the gluon-gluon production of (empty triangle) and 
B* (solid triangle) in comparison with the predictions of the fragmentation model 
(frag.) and the results of ref.[8] (curves) and ref.[10] (stars). 

Fig. 2. da/dpT at = 100 GeV for B* (solid histogram) and B^ (dashed histogram) 
in comparison with the predictions of the fragmentation model (solid and dashed 
curves for B* and B^ respectively), and the result of ref.[8] for B^. (dots). 

Fig. 3. da/dz at = 100 GeV for B* (solid histogram) and Be (dashed histogram) in 
comparison with the predictions of the fragmentation model (curves) . 

Fig. 4. da/dpT in the hadronic production of B* (solid histogram) and B^ (dashed his- 
togram) in comparison with the predictions of the fragmentation model (curves) at 
the pp-coUision energy ^/s —1.8 TeV. 

Fig. 5. The same as in Fig. 4 with the \/§ > 60 GeV cut off. 

Fig. 6. The total cross-section summed over spins of the P-levels in the gluon-gluon pro- 
duction (dots) and its fit (solid line curve) in comparison with the prediction of the 
fragmentation model (dashed line curve). 

Fig. 7. da/dz at a/^ = 100 GeV for ^Pi-, ^Pq-, ^Pi-, and ^P2-lcvcls (a, b, c and d figures, 
correspondingly) in comparison with the predictions of the fragmentation model 
(curves). 

Fig. 8. da/dpT with the same notations as in Fig. 7. 

Fig. 9. The differential cross-section of the hadronic pp production of the P-level Be mesons 
at ^/s = 1.8 TeV (histograms) in comparison with the fragmentation model results: 
^P2- solid curve, ^Pq- dotted curve, ^Pi- dash-dotted curve, and ^Pi- dashed curve. 
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